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ABSTRACT 


A rectangular cross section channel, with both a straight 
and a curved portion, was designed and built to create the 
laminar secondary Taylor-Goertler vortex flow. An aerosol 
was used to visualize these flow patterns, and photographs 
of the results are presented. A hot wire anemometer inves- 
tigation of one location in the curved section was conducted 
to obtain the mean velocity profiles and turbulence levels. 
The pressure drop along both the straight and the curved 
sections was measured and compared. A constant temperature 
wall heater was also designed and installed in the straight 
section of the channel. Approximate values of the heating 
losses were obtained and the overall and local heat transfer 
coefficients were calculated for the straight portion of the 
channel. It was concluded that the vortices develop with 
both velocity increases and distance along the channel, and 
Emde the nee aa pressure drop due to the vortices 
varied from near zero to approximately 20% with increasing 


Reynolds numbers. 


- 
a _ 
- 7 
- oe _ 
_ 
a . 7 
7 


—_ 


m >t 7 7. i 


- 2 hate! oh. 2e 


> 
ae 


mr ee 


=, ren 
Sao 

BD LP ©) Pe aerenn 
Ca aes aeee 
a aw We - 
ne se 


a th 4 &. >in 


a 


1a @ >” 


aie - um * 
| a ee); a 
» Ste ae 


2) se ee = OD Gyr 


gee ies 


» 741) a Ss 
, 1s 7 
7 o 
] 
7 
: -" 
a 

* oe Ly 

Ss: 7 
7 

7 7 


— 





TABLE OF CONTENTS 


ihe INTRODUCTION -----------~------~--~------------------- idk 
A. DESCRIPTION OF VORTICES --------------------- iat 
B. BRIEF HISTORY ~-~---~------~--------------------- iL 
ia NATURE OF THE PROBLEM ---------------~-~--~--------- Ig, 
A. INTENT OF THIS STUDY ------------------------ is 
B. DESIGN REQUIREMFNTS ~------------------------- 18 
C. DESCRIPTION OF APPARATUS ~----------~---------- 24 
III. FLOW VISUALIZATION ------------------------------ 39 
A. EXPERIMENTAL CONSIDFRATIONS ----------------- 39 
B. PATTERNS AND PICTURES OF THF FLOW ----------- 4l 
ie VELOCITY AND TURBULENCE MFASUREMENTS ~------------ 51 
A. EXPERIMENTAL PROCEDURES --------------------- pall 
B. VELOCITY PROFILES --------------------------- Ss 
C. TURBULENCE LEVELS. —------~----------- eee 62 
Wie PRESSURE DROP MEASUREMENTS ~---------------------- Va 
A. EXPFRIMENTAL PROCEDURES ~--------------------- Ve 
B. PRESSURE DROP ------------------------------- 74 
Wide, HEAT TRANSFER MFASUREMFNTS ---------------------- 85 
A. EXPERIMENTAL PROCEDURES -----~---------------- 85 
B. RESULTS ------------~------------------------- 89 
VII. CONCLUSIONS ------------------------------------- 1@i2 
VIII. RECOMMENDATIONS ----~--~---~--~------------ 104 
APPENDIX A: ERROR ANALYSIS ~--------------------------- 106 


Sebo von y Be Se COURUIER RESULTS FOR VELOCITY 


AND TURBULENCE --------------------------- 108 





PePENDIX C: CALCULATIONS AND COMPUTER RESULTS 


FOR PRESSURE DROP ------------------------ 119 
APPENDIX D: COMPUTER RESULTS FOR HEAT TRANSFER -~------ 129 
LIST OF REFERENCES ------------3------------------------ 158 
INITIAL DISTRIBUTION LIST ------------~----------------- 161 
PORM DD 1473 -<-----==--=-----------~- ~~~ =~ - - - ----- == =~ - 162 





Table 


i. 


ot. 


theo r LABLES 


PRmooUkE DROP PER FOOT FOR TWO FOOT 
INTERVALS CALCULATED BY LEAST SQUARES --------- 78 


CRs UE DROP PRR FOOT FOR ONE. FOOT 
INTERVALS CALCULATED BY LEAST SQUARES --------- 80 


PReoocme DROP PER@EOOT FOR TWO FOOT 
INC oe oes) Yao = a 81 





Figure 
1 


De. 


HO. 
LIGA 
ee 
Jb Sas 
14. 
joe 
iG. 
Lae 
aS. 
He 
ZO. 
jad We 
OL. 


Zea 


biol Or FIGURES 


Sketch of Taylor-Goertler Vortices ------------ 1D 
Reynolds Number vs. Dean Number --------------- Zi 
Cross Section of Channel -------~------~--------- 25 
Arrangement of Apparatus ---------------------- 27 
Aerosol Generating Equipment ------------------ 29 
Lighting and Camera Arrangement --------------- 31 
Hot Wire Anemometer Arrangement --------------- 22 
Pressure Drop Measuring Equipment ------------- 34 
Heating Apparatus ----------------------------- 3) fi 
Photographs of Vortices Re = 609 ------------~ 45 
Photographs of Vortices Re = 875 ------------- 46 
Photographs of Vortices Re = various --------- 47 
moO gmamns Of Vortices “Re — 10335 - [-——----—a- 48 
Photographs of Vortices Re = 1166 ------------ 49 
Photographs of Vortices Re = 1457 ------------ oe 
Ve loewty Profiles Re = 609 ------------------- 55 
Velocity Profiles Re = 875 -----------<------- 56 
Velocity Profiles Re = 1166 ------------------ 57 
Velocity Profiles Re = 1298 ------------------ 58 
Velocity Profiles Re = 1457 ------------------ 59 
Velocity Profiles Re = 1722 ------------------ 60 
Velocity Profiles Re = 1934 ------------------ 61 
Typical Position of Hot Wire ------------------ ao 





24. 
Zoe 
ZO. 
LIE 
Zo. 
eo). 
Bo. 


Bi. 


NB 


a. 


34. 
S26 
BiG. 
ae) 


Sy 3)E 


Turbulence Intensity Re 609 ---------------- 64 
Turbulence Intensity Re 875 ---r rrr cme cnn 65 
Turbulence Intensity Re 1166 --------------- 66 
Turbulence Intensity Re 2 See 2 82226 ee 67 
Turbulence Intensity Re 1457 --------------- 69 
Turbulence Intensity Re L722 -----<--<-------- 70 
Turbulence Intensity Re 1934 --------------- gat 
Pressure Below Atmospheric vs 

eee = ee Ge ee SSeS 75 


Percentage Increase in Pressure Drop 


vs. Reynolds Number -----------<-<----------- 82 
Pressure Drop Coefficient vs. 

Pressure TapS ------<-- 3-9-9930 orn enn 84 
Heat Transfer Coefficient Tw = 125 ----------- os 
Heat Transfer Coefficient Tw = 133 ----------- 96 
Heat Transfer Coefficient Tw = 143 ----------- 97 
Heat Transfer Coefficient Tw = 139 ---------=- 98 
Heat Transfer Coefficient Tw = 131] ----------- 99 





Aw 


As 


Cp 
Cp 
De 


Dh 


Go 


Hm 


Ho 


TABLE OF SYMBOLS 


area 
cross sectional area 

area of the wall heater 

area of a single heater strip 

half the height of the channel 

half the width of the channel 
pressure drop coefficient 

specific heat 

Dean number 

hydraulic diameter 

spacing of channel wall 

Goertler number 

Grashoff number 

local heat transfer coefficient 

mean heat transfer coefficient 
overall Hone ELanoret coer ricreme 
thermal conductivity 

a Aeance along the channel in the flow direction 
Nusselt number 

Prandtl number 

pressure, a relative measured value 
heat energy transferred by conduction and convection 
heat energy transferred by radiation 


Reynolds number based on hydraulic diameter 





Ri 


Ro 


Ta 
Tb 


Tm 


magic Oteiimnem wall Onecylinder 
radius of concave wall 
temperature 

Taylor number 

bulk temperature of a fluid 

log mean temperature difference 
average temperature of the wall heater 
mean velocity of flow 

volumetric flow rate 

boundary layer thickness 
emissivity 

mass density 

Stefan-Boltzmann Constant 
viscosity 


kinematic viscosity 





ACKNOWLEDGEMENT 


The author wishes to express his appreciation to all 
those who made this work possible. It is impossible to give 
a complete list of people to whom the author is indebted but 
it is hoped that by naming those who deserve a special note 
of thanks, all will know that their efforts were appreciated. 
The guidance and advice of the author's thesis advisor, 
Professor M. Kelleher, and his efforts made this undertaking 
possible. The personnel of the Mechanical Engineering Shop 
and Peet culatie Messrs. K. Mothersell, G. Bixler, and 
T. Christian deserve a special thanks for their timely and 
careful assistance in preparation of the experimental appa- 
Matus. The author 1s very grateful for the assistance of 
Mrs. V. Culley of the Mechanical Engineering Office in obtaining 
material and completing administrative requirements. Final- 
ly, to my wife Catherine, for her patience through repeated 
Manuscript typings and for her support throughout the study, a 


heartfelt thank you. 


10 


a 

ces 2 . 
7 4 
fo a | 
oa 

ee sy . 





a 


a 


tes wa 
= Vtigara: et Fe 
Sn bead eee 


7 
= 


. 


Eo 


We ae 
oe, tee 
os ee m uel 
TRG a= 2 | 
iy ate 6 a i ae 
Simson et woe 2 





et en en 


: 
: 





-_ 1497 jb Pail BG Tt 
: : 7 : ; 
Stele yo = ) 
a 
i. : a 
- - 7 a 7 
| - }: 2 7 i 
- 7. ae ie 
7 7 - 
| . nan 
. 7 7 : 
Gg -_ 5 
en es 
2 - 4 7 
7 : 7 o. 
- 


ie) ROD UC Thon 


PaeeooeoCRIPTION OF VORTICES 

The phenomenon best known as Taylor-Goertler vortices is 
@eateed by a body force instability in fluid flows. Speci- 
fically, Taylor-Goertler vortices are a secondary flow 
induced by centrifugal forces. In a channel, such as used 
in this investigation, the fluid near the center eae 
channel is subjected to larger centrifugal forces than the 
Slower moving fluid near the concave wall. Thus, the ten- 
dency is for the fluid in the center of the channel to move 
outward toward the wall. The fluid near the wall is unable 
to resist this action and must move in the spanwise direction 
and then radially inward replacing the central fluid. Once 
in the center of the channel, this fluid obtains the higher 
flow velocity, and the tendency to move toward the concave 
eee) This cyclic motion forms the rotating Tay lor-Goertler 
vortices. 

Sucihievyortlces, Caused by centrifugal femces meeceur in 
Many cases of curved flow. Taylor-Goertler vortices are 
laminar vortices with their axes in the direction of the 
main flow, and with secondary velocities in both perpendic- 
Meat Girections. Such spiral vortices occur in counter 
rotating pairs. There is frequently a regularly spaced 
cellular structure associated with these vortices. A sketch 


of this type of flow pattern is shown in Figure l. 
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Figure 1. Sketch of Taylor-Goertler Vortices. 
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Pee oRlEr HISTORY 

ine finste person to consider the instability of curved 
flows was Lord Rayleigh who found a stability criterion for 
maascid Liuid rotating symmetrically about an axis [Ref. 1]. 
His Criterion for a stable flow required that the product of 
the local circumferential velocity and the local radius of 
curvature either increase or at least remain constant as the 
radius increases. An example of an unstable flow described 
by Lord Rayleigh is the flow between a fixed outer cylinder 
peel CONnCentrically rotating inner cylinder. In~1923, 
G. I. Taylor published an extensive ae eac and experimental 
study of such unstable Couette flows of viscous fluids 
between rotating cylinders [Ref. 2]. The Taylor number, 


defined as 





_ uU*d a 
ise A 5 


was shown to be a characteristic parameter of such flows with 
eeerrtical value for the onset of the instabilaty. The 
secondary flow observed between the cylinders when the outer 
was fixed and the inner rotated has been shown to be vortices 
whose axes were in the circumferential direction. These 
vortices were named as Taylor vortices. In recent years, 
there have been a large number of investigations of cylindri- 
cal Couette flow. The torque required to rotate a cylinder 

in the presence of Taylor vortices was measured by P. Castle, 
et al. [Ref. 3]. The same type of instability in the boundary 


layer along concave walls was first investigated by 


alee: 





H. Goertler in 1940 [Ref. 4]. The Goertler number, which 
is Similar to the Taylor number, is a parameter characteris- 


tic of curved boundary layer flows and is defined as 





* 
Go = ¥ : a 
Vv R 
where 6 is the boundary layer thickness. When the Goertler 


number exceeds a critical value, the presence of Goertler 
vortices is observed. H. W. Liepmann [Ref. 5] was among the 
investigators to continue and verify Goertler's experimental 
work. The approximate analytic results obtained by Goertler 
have been verified with an exact solution by the German 
investigator G. Hammerlin as reported by H. Schlichting 
[Ref. 6]. These results were also verified with an extensive 
Memerical solution completed by A. M. O. Smith [Ref. 7]. 

The problem of instability of viscous fluids flowing in 
a curved channel was first considered by W. R. Dean in 1928 
[Ref. 8]. Dean found a parameter, similar to the Taylor and 
Goertler numbers, which governs the onset of Taylor-Goertler 
vortices in ie cAteeRe: The analytic work done by Dean 
has been verified by W. H. Reid [Ref. 9] among others. The 
experimental work on secondary flows in channels seems very 
limited. Until just recently, the possible effects of 
Taylor-Goertler vortices on heat transfer had not been 
studied. In 1965, L. Persen [{Ref. 10] considered the effect 
of Goertler type vortices on the heat transfer froma wall 
for the special cases of both a very high and a very low 


Prandtl number. The first experimental work with the effect 
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of Taylor-Goertler vortices on heat transfer in a boundary 
layer was reported by P. McCormack et al. [Ref. 11] in their 
1970 study. With some limiting assumptions, a numerical 
solution for the heat transfer in curved rectangular chan- 
nels was reported by K. Cheng and M. Akiyama in 1969 

[Ref. 12]. Analytic and experimental results for a fully 
developed, constant wall heat flux, square cross section 
curved channel were obtained in 1970 by Y. Mori, et al. 
Peet. 13). 

Although Taylor-Goertler vortices are a laminar flow, 
they have secondary velocity components and are more complex 
than one would expect. Furthermore, as indicated by H. W. 
Liepmenn [Ref. 5], these vortex flows are potentially a key 
to understanding transition to turbulent flow. Taylor- 
Goertler vortices have similarities to other vortex flows 
such as wing tip vortices and the vortex rolls in forced 
convection heating of fluid layers described in the paper by 
Ms Akiyama et al. [Ref. 14]. The striations seen at stag- 
feiton points on Dee bodies and the cross hatching on 
reentry bodies have also been attributed, at least in part, 
to the presence of Taylor-Goertler vortices, as reported by 
MeerPersen [Ref. 15]. 

There are many possible applications that could result 
from investigation and understanding of these vortex phenom- 
enon. Potential application include improved cooling of 
turbine blades and other external surfaces, reduction of 


pressure losses in bends, and improved mixing in laminar 


iB 













ae 
oa 7 oa _ : 
oe. o os < so 


et, 
















flows. Heat exchangers which could utilize the low pressure 
drop in laminar flow as well as the enhanced heat transfer 


in the presence of vortices would be a significant advance- 


Mente in heat transfer. 
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ii) NEC, Or SLE PROBLEM 


ieeeLNIENT OF THIS STUDY 

Because future application of the Taylor-Goertler vortex 
phenomenon may take place in closed channels, it was intended 
to study the effect of these secondary flows in a curved 
channel of rectangular cross sectional area. To investigate 
Taylor-Goertler vortices in a curved rectangular channel, the 
presence of the vortices had to be verified. It was intended 
to visualize the flow in order to confirm the presence of 
the vortices. How rapidly the vortices would form and to 
what extent they would fill the channel was not known. It 
was intended to use the flow visualization technique to 
study, as far as possible, the size, the form, and the 
development of the vortices. 

Other planned measurements included a hot wire anemometer 
check of the turbulence level to insure that the flow was 
Taminar. The Pe oaiey profiles in the channel were also to 
be obtained to determine what effect the secondary flow would 
have on the mean velocity. The additional velocity com- 
ponents in Taylor-Goertler vortices produce Reynolds stress 
not present in other laminar flows. It was planned to 
measure the pressure drop along both the straight and 
curved portions of the channel used in this investigation 
to determine the effect, if any, of these Reynolds stresses 
on the pressure drop. An increase in the pressure drop 


would be anticipated because it is known that additional 
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energy is required to maintain the secondary flow. This 
energy 1S supplied to the Taylor-Goertler vortices by an 
interaction of the mean and secondary velocity components. 
Mm@ew works Of G. I. Taylor [Ref. 2] and P. Castle et al. 
[Ref. 3] have shown the increase in torque required to 
rotate a cylinder in the presence of Taylor vortices. 

The Taylor-Goertler vortices would also be expected to 
effect the heat transfer from the concave wall. This effect 
would result from the secondary velocity components trans- 
porting heated fluid from the wall towards the center of the 
channel, and carrying cooler fluid from the center to the 
wall. It was intended to study the requirements for measur- 
ing the effects of this action on the heat transfer from 


the concave wall. 


Pee oeolLGN REOULREMENTS 

In order to study Taylor-Goertler vortices, an apparatus 
had to be built which would give the values of the non- 
dimensional paraméters required to produce the secondary 
Flow. The limit of stability of the basic laminar flow in 
a curved channel against the formation of Taylor-Goertler 
vortices is well known, and is presented as a critical Dean 


number by W. Reid [{Ref. 9]. If the Dean number, defined as 


_ u*a d 
De a \) Rp i 





exceeds a value of 36, it is presumed that Taylor-Goertler 


vortices will be present. Furthermore, it is known that the 
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rate of development of the vortices will be somewhat depen- 
dent on how much greater than 36 the Dean number becomes. 
It is known that the size of the vortices will depend in 
some manner on the dimensions of the channel. Investigators 
of Taylor-Goertler vortices on concave walls have reported 
vortices from the size of the boundary layer to ten times 
larger, dependent on geometric factors [Ref. 7]. The Dean 
number had to be greater than the critical value of 36 and 
be adjusted within a workable range. One hundred fifty was 
taken as the upper limit of this range, as it was assumed 
the vortices would be fully developed by at least this value. 
The minimum critical Reynolds number, based on hydraulic 
diameter, at which turbulent flow is observed in a parallel- 
plate channel is given by G. Beavers, et al. [Ref. 16] as 
Re = 2200. Several investigators including H. Liepmann 
[Ref. 5] and R. Nunge [Ref. 17] have reported that the 
decrease in the critical Reynolds number due to curvature is 
small, especially for a large radius of curvature. There- 
fore, the maximum Reynolds number, based on hydraulic dia- 
meter, selected for this investigation was Re = 2000. The 
Reynolds number had to be high enough for the flow velocities 
to be measured accurately. A minimum of Re = 100 was used 
for the design of the experimental apparatus. 
The design procedure involved checking the effects of 
several factors on the range of the nondimensional variables. 
The working fluid, the channel height, and the flow rate 


selected, each have an effect on both the Reynolds number 


iS 





and the Dean number. The radius of curvature also effects 
the Dean number. Trial and error was used to select the 
height and width of the channel's cross section, and the 
following inequality in the Reynolds number was checked to 


give the range of velocities 
NOG = Re < 2000 


The inequality in the Dean number given below was then 


checked to give the range of the radius of curvature. 
Some = ES 0 


When the velocity range was practical to obtain and easy to 
measure, and when the radius of curvature was reasonable for 
construction, the design variables were considered feasible. 
The resulting relationship between the Reynolds number, 
based on hydraulic diameter, and the Dean number with the 
mietoror radius of curvature to plate spacing as a parameter 
1s shown in Figure 2. The selected design is represented by 
the line R/d = 48. 

The design of the apparatus for this investigation was 
also constrained by time and cost limitations. The material 
and equipment to be used had to be workable, available with- 
out long delay, and inexpensive. Water and air were the 
only working fluids considered as practical. The flow 
velocities at which water satisfied the Reynolds number and 
the Dean number requirements were too low to measure by 
ordinary means. Furthermore, water requires more elaborate 


construction to prevent leakage from the channel. For these 
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reasons, air was selected as the working fluid. For design 
purposes, the value of the fluid properties at estimated 
operating temperatures were taken from reference tables. 
The length of the entrance region in the design channel was 
checked to see that it would not be excessive. To insure 
smooth entering flow, the entrance contraction section was 
designed in accordance with considerations given by M. Cohen 
ama. Ritchie [Ref. 18]. The intent to visualize the flow 
required the use of a workable, transparent, inexpensive 
material for the walls of the channel. Plexiglas was 
selected for this purpose. 

It was hoped that the results obtained for both pressure 
drop and heat transfer would be comparable to infinite 
parallel plate solutions, so a large aspect ratio was chosen. 
The pressure drop for flow between parallel plates is known 
meomedn exact analytic solution and is not of a large magni- 
tude. This magnitude was compared with the sensitivity of 
‘the laboratory micromanometer available. It was decided 
that the pressure drop could be measured, but only if a 
sufficiently long test section was built. No exact solution 
exists for the pressure drop in a curved channel in the 
presence of Taylor-Goertler vortices. However, the studies 
by R. Nunge [Ref. 17] and by K. Cheng and M. Akiyama [Ref. 12] 
indicates a rather small increase in the friction factor in 
curved channels of large aspect ratio. It was assumed that 
the curved test section would be sufficiently long to allow 


measurement of the pressure drop. 
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To measure the effects of Tee oe vortices on 
the heat transfer from a constant temperature concave wall 
of a curved channel requires comparative data for one wall 
Srmoaecilat plate channel. Heat transfer data for a flat 
plate channel with one constant temperature wall is very 
limited. Limiting values of the Nusselt number and other 
comparative values can be obtained from R. Shah and A. Lon- 
don's report [Ref. 19]. The design of a heater plate to 
measure this type of heat transfer presented special require- 
ments. As with any flat plate heater, more power would be 
convected away from the upstream edge than from points 
farther downstream. Therefore, to obtain an uniform wall 
temperature, by electrical heating, separate resistance 
elements with individual voltage controls were reguired. A 
multiple channel regulated power controller was not immedi- 
ately available, but it was determined that such an instru- 
ment could be constructed. Grade A nichrome was selected 
for the heating elements because of this metal's large 
specific resistance, small temperature coefficient of resis- 
tance, and relative on cost. The many power leads and 
temperature recording leads would have to reach the heating 
element through the channel wall. For this reason a castable 
plastic resin was selected for constructing the heater sec- 
men. Although only a flat heater plate was to be constructed 
for the first tests, design considerations were given to 
adapting this construction method to a heater plate for the 
curved section. Other aspects of the heated section con- 


Peat lon wllilebe discussed Later. To insure that the heat 
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maeeanscter Gata would be for forced convection, and would not 
@Gemtaln the effects of free convection, the value of the 
governing dimensionless group was checked against the 


requirement 


Gr * Pr * a < 1000 


Ganyven by J. Holman [Ref. 20]. It was found that this 
requirement imposed a significant temperature difference 
limitation for water as the working fluid, but allowed a 
large temperature range for the selected working fluid, air. 
The effects of the many design requirements discussed 
above was to make the design problem somewhat difficult. 
However, the working fluid, the geometry, and the flow 
velocities to be studied are of values frequently encountered 


PimeGactice. 


See SoCAIPTION OF APPARATUS 

The cross section of the channel was 0.250 inches high 
and 10.0 inches wide. It had an aspect ratio of 40 and a 
cross sectional area of 0.017361 square feet. The hydraulic 
diameter of this channel was 0.040650 feet. The channel was 
made of two one-quarter inch thick sheets of plexiglas 
separated by rolled metal spacers which also formed the 
sides of the channel as shown in Figure 3. Steel was used 
for the spacers in the straight section of the channel and 
for added flexibility, aluminum was used in the curved por- 


muon. A Straight portion of this channel, consisting of a 


24 





“TOUURUD FO UOTISeS SSOID “°fE SanbTY 


C2 PLTPOS 


3zoddns poom * 


ESSSSSS SSSNSS | 
RRL LLL LILI ILL L IL ELLE LIL LLL LLL LLL 

: (2 =e ves See ee eee eee 

bee ae a TE LO EA PDD 


‘ i 
SL" .0°0T uae 


zeoeds Tezow seT6txeTd 





two foot long entrance section and a two foot long test 
section, preceded the curved section. The radius of curva- 
ture of the interior concave wall of the curved section was 
one foot. The curved section was formed by heating the 
plexiglas and bending it over a carefully constructed wood 
frame. This wood frame was also used to support the entire 
channel. The curved portion of the channel completed a 
180 degree turn and was followed by a short straight section. 
The flow of air entered the channel through a contrac- 
tion section attached to the beginning of the entrance region 
and covered by a cheesecloth screen for most of the experi- 
mental work. The flow left the channel through an exhaust 
nozzel which was connected to flexible tubing of one inch 
inside diameter. The tubing led the flow through a Fischer 
and Porter Company variable area flow meter, model number 
10 A3565. This rotometer had a 100% full scale flow rate 
of 11.1 standard cubic feet of air per minute and an accuracy 
of +0.5% of full scale. The flow was drawn through the flow 
meter, the lemibile tubing, and the channel by an electri- 
cally driven Cadillac, model G12, centrifugal blower. The 
blower Ba, ds controlled by varying motor voltage with a 
variac. The voltage was supplied by a Sorenson, R 1050, 
A.C. voltage regulator. The general arrangement of this 
apparatus is shown in Figure 4. 
The idea of using an aerosol for visualizing the flow 
was taken from a paper by O. Griffin and C. Votaw [Ref. 21]. 


The advantages and ease of using an aerosol are derived from 


Z6 








uoT OES 
4FS04 
peaano 


a a 


rg ae ial! 





il 


*>snqzeireddy 


uOoTIOSS 


ic 


3S03 3UbTAAS 


JO qZuewshbuezzy «°P Sanbt J 


aTzzou ysneyxs 





Sea eiiileweo: 


ujZbueT souer US 


DeTICA ee 
UO) 


IOMOTO 
uot zons 


uoTzaANS 
uoTzOeAAUOD 


27 





the non-toxic and non-corrosive properties of the DOP, 

di (2-ethylhexyl)-phthalate, chemical used. The equipment 
used to make the aerosol was built in accordance with a 
United States Naval Research Laboratory report by W. Echols 
and J. Young [Ref. 22] and personal communications with 

O. M. Griffin. Air from a regulated and filtered supply was 
piped to an atomizer nozzle submerged ina jar of DOP. In 
addition to the pressure regulator, there was a throttling 
valve on the air line and the depth of submergence of the 
nozzle could be adjusted. The aerosol from the top of this 
jar was collected and led to a jet impactor in a second jar. 
The effect of the jet impactor was to remove excessively 
large particles and to make a smaller, more uniform aerosol 
particle size. The aerosol was then led to a settiing 
chamber where the pressure was reduced to nearly atmospheric 
level by bleeding off the excess aerosol. Flexible injec- 
tion tubing of 0.250 inch inside diameter connected the 
settling chamber to the contraction section where the 
aerosol entered the channel. The aerosol generating equip- 
rove pictured in Figure 5. The character of the aerosol 
used was such that visualization of the flow was facilitated 
by illuminating slits along the channel and viewing them at 
approximately 90° to the light. A single 650 watt Colortran 
light, model 100-071, was placed on the side of the channel 
opposite the viewer. The light was collimated with a paper 
cone lined with aluminum foil. Flat black paper was taped 


to the outside of the channel walls to make four narrow 
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Figure 5. Aerosol Generating Equipment. 
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slits. The light beam was directed through one of these 
slits so that a small area of the flow could be visualized. 
A Nikomat 35 mm camera and a Nikon 55 mm £3.5 Macro-Nikkor 
lens was used to obtain photographs of the flow patterns. 
A sketch of the lighting and camera arrangement is shown in 
Figure 6. The visualization techniques and procedures will 
be discussed more fully when the photographs are presented. 
In addition to visualizing the flow it was planned to 
measure the velocity profile. Due to the size of the chan- 
nel, a sub miniature hot wire anemometer probe, Thermo- 
Systems Inc. model 1279, which had a 90° bend in the supports 
was used to measure the turbulence levels. This probe was 
located three inches from the center line of the channel and 
20 inches of arc length downstream from the start of the 
curved portion of the channel. The probe was placed ina 
plug and inserted through the convex wall of the channel. 
The inside surface of this plug fit flush with -the wall of 
the channel, and had the same radius of curvature. The 
02025 inch STO a body of the hot wire anemometer protruded 
from the plug and was fastened to a micrometer barrel which 
was mounted on a machined block attached to the outside of 
the channel. In this way, the movement of the wire from 
very close to the concave wall of the channel to near the 
convex wall could be Boe ene and measured. This arrange- 
ment is shown in Figure 7. The leads from this probe were 
connected to a Thermo-System Inc. hot wire anemometer 


system, model 1050, including a Hewlett-Packard RMS Meter. 
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Figure 7. Hot Wire Anemometer Arrangement. 
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A Thermo-Systems Inc. calibrator, model 1125, was used to 
calibrate the hot wire anemometer. 

In order to measure the pressure drop along the length 
of the channel, pressure taps and a micromanometer were 
required. Pressure taps were located every three inches 
along the entire length of the channel, starting two inches 
from the beginning of the entrance section and alternately 
set 0.250 inches to the right or left of the center line. 
Along the curved section, the three inch spacing was 
measured on an arc of one foot radius, resulting in a total 
of 28 pressure taps with number 17 being two inches into 
the curved section. These pressure taps consisted of a 
0.040 inch diameter hole drilled into the channel with a 
Sa@imeplexiglasjtube glued into a 0.250 inch hole counter- 
bored over the smaller holes. Flexible tubing with pressure 
fittings was connected to these short tubes to facilitate 
connecting and disconnecting the micromanometer. The instru- 
ment used in this investigation was an E. Vernon Hill and 
Company Type C ie comenenevee with a range up to two inches 
of water and a sensitivity of 0.001 inches of water. This 
micromanometer is of the re-zeroing type where the reservoir 
is raised to return the meniscus to the same level it had 
before the pressure was applied, and the amount the reser- 
voir is raised indicates the pressure. A photograph of the 
channel and the pressure measuring apparatus is shown in 


Figure 8. 
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The bulk temperature of the air was measured six inches 
after the entrance of the channel, and three inches after 
the heating plate, by Pours ih emiacete We Spaced across the 
channel at each location. Copper-constantan thremocouples 
were used. They were located at different heights in the 
channel and connected in parallel to give the average bulk 
temperature. An ice bath reference junction was provided, 
and these thermocouples and the others used in this experi- 
ment were connected through a switch to a Dymec digital 
voltmeter, model 2401B. These thermocouples and a sample 
of the others used in this study were calibrated in a 
Rosemont Laboratory Calibration System. In this constant 
temperature bath, the temperature was determined by a 
platinum resistance thermometer. 

A wall heater was constructed for the straight test 
section of the channel. This heater was formed of forty 
Mmeoeeinch thick, 0.125 inch wide strips of nichrome run 
across its 11.5 inch width. These metal ribbons were 
placed in a mold by hand and spaced approximately 0.010 
inches apart. The streamwise length of the plate was 
00 inches. A copper-constantan thermocouple was welded 
to the center of each of the nichrome strips. Power leads 
were connected to one side of each strip four inches from 
pre center of the plate. On the other side of the center- 
line, the strips were connected together in pairs to make 
20 electric circuits. The 20 heaters were individually 


Supplied by a specially constructed, regulated 
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high current D.C. power supply and a 20 channel power con- 
troller. The D.C. power supply had a 200 amp capacity and 

a maximum RMS ripple of one millivolt. The power controller 
was a 20 channel potentiometer biased series regulated 
voltage source with a temperature compensated reference 
voltage. The components of this power controller were 
mounted on a common heat sink in order to improve the 
stability of operation. The power supply and controller 

are shown in Operation in Figure 9. When all the power 
leads and thermocouple leads were in place, APCO 210 epoxy 
resin with APCO 180 hardener was poured into the mold. This 
thermal set plastic was cured at room temperature for 
approximately 24 hours. When the heater plate was removed 
from the mold, it was found to be warped. A metal frame was 
constructed, securely attached to the outer edges of the 
e250 anch thick casting and adjusted to produce a Monae coals 
flat plate. A further problem experienced with the casting 
involved some of the nichrome strips near the trailing edge 
of the plate that had pulled away from or sunk into the 
plastic during drying. This resulted in rough spots on the 
surface of the plate. The loose strips were cut and removed 
from the heater, and the low spots, from the 33rd strip to 
the end of the plate, were filled with epoxy. The plate was 
then lightly sanded, resulting in a smooth surface. The 
first 16 heating circuits remained useable, which resulted 
in a heated area 8.125 inches wide by 4.384 inches in the 


Mmiow direction. 
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After the other experimental work was completed, a 
5.50 inch piece of the upper wall of the channel, between 
the two foot ten inch station and the three foot 3.50 inch 
station, was removed. The heater was placed into this cut, 
becoming part of the straight wall of the channel with the 
nichrome ribbons exposed to one flow of the air. The top 
of the heater and the surrounding channel walls were then 
covered with a 0.125 inch thick layer of Johns~-Manville 
Company Min-K insulation. The Min-K was covered with a one 
inch layer of glass wool. Each piece of the apparatus was 
cleaned and carefully prepared before the experimental work 


began. 
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Crt OW ViIovAbL ZATION 


Ae EXPERIMENTAL CONSIDERATIONS 

The aerosol generating equipment was found to be very 
versatile as to the amount and quality of the aerosol pro- 
Guieed. sAfter much trial and error, good quality aerosol was 
found to be that of relatively esrnesiia particle size and low 
exhaust pressure. This aerosol appeared much like cigarette 
smoke when discharged into the room. A supply air pressure 
of 30 pounds per square inch, or slightly less, an atomizer 
nozzle submergence of approximately 0.50 inches, and a low 
settling chamber pressure were the conditions which produced 
this aerosol. Producing this quality aerosol and getting 
the right amount of it into the channel proved to be impor- 
tant in visualizing the flow. Viewing the aerosol at the 
proper angle relative to the incident light was also impor- 
tant. For lighting, the best arrangement was to place a 
Single, strong light near the center of curvature of the 
channel and to shine it radially outward. The collimated 
beam of light was directed through narrow slits formed from 
black paper taped to the channel. If this small area of 
the flow was viewed from an angle tangent to the channel, 
the vortex patterns could be readily seen. This viewing 
location was the same as the camera location shown in 
Figure 6. Four such narrow slits were made along the curved 


feeiam Of the channel. The first slit was located 7.50 
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inches from the beginning of the curved portion of the 
channel. The second was 12.00 inches downstream from the 
first, and the third was 11.50 inches farther downstream. 
The last slit was located at the end of the curved section, 
Bpeeeanches from the third slit. At the location of this 
fourth slit, the channel enters the final straight section. 

To obtain photographs of the flow patterns in each of 
these slits, the camera with the 55mm macro lens was mounted 
on a tripod and moved as close to the channel as possible. 
The camera was aligned with the axes of the vortices by 
placing the camera on a tangent to the channel. In the case 
of the first and second slits, the flow was coming toward 
Mmiemedamera, and in the case of the third and fourth slits, 
the flow was away from the camera. Due to the many reflec- 
tions observed from the surfaces of the plexiglas, the room 
was darkened while the photographs were taken. A Gossen 
“Luna-Pro" exposure meter with a variable angle spot meter 
attachment was used to assist in obtaining the proper 
exposure sotbncse Under these difficult conditions, several 
attempts were required before satisfactory photographs were 
obtained. | 

Other arrangements for viewing the flow were also tried. 
Looking radially inward through the channel toward a less 
intense light located near the center of curvature showed a 
emetterent type of pattern. The patterns seen were streaks 
eeeiines in the direction of the flow. These patterns were 


similar to those seen by other investigators including 
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L. Persen [Ref. 25]. Several methods of injecting the aero- 
sol into the flcw were also tried. The methods involving 
small tubes or fittings did not work because sufficient 
aerosol would not pass through these passages. The results 
obtained from the other methods were for the most part the 
Same as for the method chosen. In a few cases, the early 
development patterns appeared differently, depending on how 
the aerosol entered the channel. The fully developed vortex 


pattern, however, always appeared the same. 


Pee f©FALTERNS AND PICTURES OF THE FLOW 

The pictures of the vortex patterns are presented in 
Figures ten through 15. Below each picture is a slit number, 
a Reynolds number and a Dean number indicating the location 
of the slits, the flow rate and the value of the character- 
istic parameter for secondary flow respectively. For the 
channel used in this investigation, the Dean number is 


Gebrated to the Reynolds number by 
De = 020729735 % "Re, 


The line below the vortex patterns in these enlarged pictures 
is the concave Wall of the channel. The line of reflected 
imgeiit above the Wortices is the inner wall of the channel 
which is 0.25 inches from the concave wall. The photographs 
taken at slit four appear slightly different than the others 
because the lighting angle and therefore the reflections 

from the plexiglas were slightly different. The pictures in 


each figure are arranged with slit one at the top of the page 
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and the downstream slits in order below the first. The 
Reynolds numbers of the pictures in each figure are intended 
to be of the same magnitude. There is not complete repre- 
sentation of the same Reynolds numbers or an even sampling 
of the slit location because of the moderately high failure 
meeenwot Obtaining pictures. In spite of this difficulty, 
the photographs presented are indicative of the Taylor- 
Goertler vortices as observed. The figures are in order of 
increasing average Reynolds number. 

A number of distinct patterns appear repeatedly in these 
pictures. At the lowest velocities, the aerosol remains 
near the concave wall in a pattern that is not clearly a 
fully developed vortex flow. Examination of these patterns 
in Figures ida, 10b, lla, and 12a show that the ends are 
curled up. This same pattern has been observed farther 
downstream in slit three at flow rates lower than pictured 
here. This pattern occurs predominately near the beginning 
of the curved section. From extended observation of both 
this pattern and the motion at the curled end, it is felt 
meee they represent the beginning of a developing weak vortex. 

The next form considered looks much like a "mushroom." 
This pattern is exemplified by Figures 12b, 11b and the right 
Side of 10c. The appearance of two different patterns under 
the same conditions, as in Figure l10c, was observed occasion- 
ally. Such occurrances are possibly explained in terms of 
the actual finite width of the channel which produces a 


stronger velocity gradient near the side walls than in the 
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middle. A difference in the placement of the aerosol injec- 
tion tubes could also contribute to the appearance of dif- 
ferent patterns at the same flow rate. Furthermore, as the 
vortices develop and occupy more space in the channel, some 
will be distorted to make room for the others. At higher 
Reynolds numbers, the top of the curl pattern, as in 

Figures l3a, 14a, and 15a, could be considered a "mushroom" 
on a small scale. The "mushroom" shape occurs farther down- 
stream or at higher velocities than the first pattern 
discussed. This shape always proceeds in distance or flow 
rate, the first clear vortex patterns. The "mushroom" has 
been observed to lead directly into distinct vortices and 
therefore should be considered as a developing vortex. 

The next pattern is the round form shown in Figures llc, 
13b, 14b, and 15b. These forms represent a counter rotating 
pair of vortices with flow leaving the concave wall in the 
center and curling back toward the wall on either side. 
These round vortex patterns do not appear to occupy all the 
Space in the eae. In Figure 13b there appears to be 
space above the vortices, and in Figure 15b the space appears 
below the vortices. However, the aerosol which is what was 
photographed does not necessarily show the entire vortex 
pattern. The aerosol shows a streakline line pattern which 
may be somewhat dependent on how the aerosol enters the 
developing vortices. As the flow rate and the downstream 
distance increases, the fully developed vortex patterns 


appear. These vortices, as shown in Figures 10d, lld, and 
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2c areythe full height of the channel and occupy a nearly 
collulapeareca. ouch patterns could be considered as fully 
developed Taylor-Goertler vortices. 

As the velocities continued to increase, the vortices 
grew in width until they began to crowd one another, as seen 
mmemeeures 13c, 4c, 14d and 15d. If the distance or the 
Reynolds number were increased farther, an unsteady oscil- 
latory motion developed in the spanwise direction. This 
action which could be observed by eye, was difficult to show 
a a Still photograph. Figure 13d is the result of such an 
attempt. The observations indicated that the sideways oscil- 
lations of the vortices may involve the absorption of one 
vortex by an adjacent growing one. It was apparent that an 
increase of oscillatory motion would have led to turbulence. 

It was evident from the photographic study that the 
vortices did require some distance to develop. Several 
different flow patterns were observed during the development. 
At high-flow rates, less distance was required for the 
vortices to become full size. There were several complex- 
Lies Sloseree in the flow patterns which could lead to 
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Figure 10 
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Figure 14 
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IV. VELOCITY AND TURBULENCE MEASUREMENTS 


Pee ee RIMENTAL PROCEDURES 

The sub miniature hot wire anemometer probe to be used 
to obtain the mean velocity profiles was connected to the 
anemometer. The cold resistance of the wire was measured 
as 7.90 ohms. Based on an overheat ratio of 1.5, the 
anemometer's operating resistance was set at a value of 
11.85 ohms. The anemometer's signal was fed into an oscil- 
liscope so that the patterns could be observed. Stability 
and trim controls were adjusted for maximum frequency 
response. The probe was centered in the chamber of the 
calibrator. A filtered air supply was connected to the 
calibrator, and the pressure drop across the calibrator's 
flow nozzle was measured by an incline manometer. With the 
anemometer switched to run, the flow rate through the 
calibrator was adjusted to different levels. The pressure 
drop in inches Be eicex, the bridge output voltage, and the 
RMS voltage were recorded. This data was reduced by the 
calibration program shown in Appendix B and the square of 
the output voltage was plotted as a straight line against 
the square root of the flow velocity in feet per second. 
The slope and intercept of that straight line were the 
values required to reduce further data, and were found to 


be 0.710 and 2.020 respectively. 
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The hot wire anemometer was removed from the calibrator 
and inserted as far into the plexiglas mounting plug as 
possible. The mounting block for the micrometer barrel was 
meeaecneda to the stem of the probe, and the distance from 
the base of the block to both the top of the plug and to 
the hot wire were measured. The difference in these machine 
vernier readings was 0.050 inches. This was as close as 
the hot wire anemometer could come to the convex wall of the 
channel. The micrometer barrel and a straight piece of 
steel wire had been used to measure the distance across the 
channel at the location of the plug. This distance was 
founda to be 0.240 inches, indicating that the probe would 
have a travel of 0.190 inches across the channel. When the 
plug was inserted into the channel and the micrometer was 
attached to the mounting block, the micrometer reading was 
0.460 inches. The hot wire anemometer was carefully moved 
until the supports for the wire contacted the concave wall 
of the channel. The micrometer reading at this point was 
0.270 inches, verifying the 0.190 inches of travel. The 
wire itself was very close to, but not touching, the wall. 
The data was taken between these micrometer settings of 
0.270 inches and 0.460 inches. 

Before the data was taken, the channel, the blower, and 
the electronic instruments were allowed a warm up period of 
about one hour. After the warm up, the flow rate was 
adjusted and a minimum of 30 minutes was allowed for steady 


state to be achieved. At each flow rate, the temperature of 
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the incoming fluid was measured and the micrometer was moved 
0.010 inch between data positions. At each position, the 
flow rate was checked to see that it was stable, and the 

Me ometer setting, the output voltage, and the RMS voltage 
were recorded. This data was reduced by a computer program 
which plotted the mean velocity and the turbulence intensity 
against the distance from the concave wall. This program 


and the printed results are contained in Appendix B. 


Pee veEGOCITY PROFILES 

The velocity recorded and shown in the figures which 
follow was not just the mean velocity component down the 
channel, but the vector sum of the streamwise and the radial 
velocity components. The portion of the response that was 
due to the radially inward or outward flow could not be 
determined by this one set of measurements. However, at the 
flow rates of this experiment, a nearly parabolic velocity 
profile would be expected if the secondary flow was not 
present [Ref. 17]. This has also been shown for curved 
channels of square cross section [Ref. 13]. The position 
of the hot wire relative to the vortices was not determined 
during the tests. Considering the nature of vortex motion, 
Peomert in the position of the vortex relative to the probe 
could have made a difference in the velocity indicated. The 
location of the hot wire anemometer was very close to the 
Second slit used in flow visualization. Thus, for each flow 
rate, the vortices may be considered to have developed to a 


corresponding level. 
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Figures 16 through 22 show the mean velocity in feet per 
second versus the distance from the concave wall in hun- 
dredths of an inch. There was a maximum uncertainty of 
approximately 0.010 inches on the location of the wall. The 
uncertainty in the velocity measurement using the T.S.I. 
system was at most 0.05 feet per second. Figure 16 shows 
a parabolic shape with a small increase in velocity indicated 
near 0.060 inches. This point may be associated with an 
active point of a developing vortex. The step in the velo- 
city shown near the wall could be associated with the probe 
actually contacting the wall, or with the error in position- 
ing. In this case, and in all other graphs, the curve should 
be extended to zero at the wall. Although a vortex was most 
probably present, Figure 17 shows no irregularities in the 
velocity profile. Perhaps the probe was located near the 
center of a vortex at this time. The same comment could 
be true of Figure 18. The magnitude of the velocity is 
increasing as would be expected. 

Figure 19 shows a variation from the parabolic profile. 
As indicated by Figure 15b, the vortices were growing at 
this Reynolds number and a different portion of the vortex 
may have been acting on the wire. In Figure 20 the velocity 
profile appears much different. In this case, the radial 
velocity component was larger near the wall than in the 
center of the channel. Figure 21 is much the same as 
Figure 20. A different profile appears in Figure 22. This 


last profile could be explained in terms of the hot wire 
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moving through the edge of a vortex. Near the wall the 
outward radial velocity would increase the mean velocity 
measured. In the center of the vortex there may be little 
effect on the hot wire anemometer. Near the top of the 
vortex the radial component of the velocity would again 
affect the probe. At the other edge of this vortex, the 
inward movement of the fluid would produce the same velocity 
profile. The length of the hot wire, 0.063 inches, would be 
sufficiently small to traverse the vortices in this manner. 
A schematic drawing of the probe in a typical position 


relative to the vortices is shown in Figure 23. 


fee LURBULENCE LEVELS 

The turbulence intensity was defined as the root-mean- 
square of the velocity fluctuation divided by the mean 
velocity. These values were calculated by the data reduc- 
tion program from the RMS meter readings and the bridge 
output voltages. The results were plotted against the 
distance from the’concave wall in Figures 24 through 30. 

The turbulence level was quite low in the center of the 
channel and higher near the walls as result of the low 
mean velocity near the walls. 

The change in the turbulence level near the wall in 
Figure 24 was due to the step change in the velocity profile 
for that location and Reynolds number. Otherwise, Figures 
mwe25, and 26 are self explanatory. In Figure 27, the 
curve representing the turbulence intensity is not as smooth 


as the other curves. Large increases in the fluctuations of 
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Figure 23. Typical Position of Hot Wire. 
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the RMS meter started at this flow rate, resulting in an 
increased uncertainty in the level of turbulence. With the 
RMS meter reading at the 0.00350 volt level, the value would 
suddenly increase as much as 0.010 and drop back again. 
The frequency of this random increase was of the order of 
30 cycles per minute. One possible cause of these large 
jumps in the local velocity would be the sideways oscillatory 
motion of the vortices. During such sideways motion, dif- 
ferent portions of the vortices could cross the hot wire 
anemometer causing the fluctuations. These fluctuations 
were present during the measurements that led to Figure 28. 
This figure represents a uniform turbulence level across the 
channel. In Figures 29 and 30, the level of the turbulence 
was larger than in the previous cases. The magnitude of 
the fluctuations was also doubled or larger. The occasion- 
ally high value of the turbulence intensity would indicate 
Euepulent flow. In other words, there was bursting to 
turbulent flow at the higher flow rates. 

The growth ana motion of the vortices within the channel 
has been indicated by the mean velocity profiles and turbu- 
lence ec ore sentad here. Evidence for the transition 


to turbulent flow has also appeared. 
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Ve PRESSURE DROP MEASUREMENTS 


ee EXPERIMENTAL PROCEDURES 

To measure the pressure drop along the channel required 
a great deal of care. Each of the 28 pressure taps were 
cleaned and checked for tightness with high pressure air. 
The flow rate in the channel was adjusted and allowed to 
Stabilize for approximately one hour before data was taken. 
The micromanometer was positioned on the table, leveled, 
and filled with the special fluid provided by the manufac- 
eater. 

The meniscus of this instrument was viewed under a 
Magnifying glass and carefully adjusted so that the curve of 
the meniscus just touched the hair line on the glass tube. 
This was established as the zero position. It was found 
that the movement of even a small weight on the laboratory 
table would effect the level and therefore the zero of this 
micromanometer. oe micromanometer was zeroed, no 
weight movement was permitted and the investigator did not 
touch the table while measurements were taken. 

The pressure fittings for each of the pressure tapS were 
numbered and taped to the edge of the channel. The pressure 
drop between pairs of pressure taps at one foot intervals 
was Measurable, but contained too much uncertainty to be of 
use. Typical values of 0.007 inches of water were measured 


as the pressure drop per foot. However, values of 0.001 
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inches higher or lower could be obtained in the same posi- 
tion of the channel at the same flow rate. It is possible 
that part of these variations were due to small fluctuations 
in the flow rate caused by the centrifical blower. Other 
causes were considered, but such fluctuations could not be 
positively identified or eliminated. 

The difference between the local pressure at each tap 
and atmospheric pressure proved to be more useful. At 
successive locations down the channel the magnitude of this 
pressure difference became larger and hence the relative 
error was less. In taking each pressure reading, the zero 
of the manometer was checked and then the flexible tubing 
from the micromanometer was connected to the pressure fit- 
ting. The pressure tap number was checked and the inlet 
temperature reading was taken and both were recorded. The 
manometer reservoir level was adjusted with the fine screw 
dial until the meniscus of the fluid had returned to its 
zero location. The rotometer was checked to insure that 
the flow rate was eran and stable. The meniscus of the 
fluid was then rechecked and if it had remained in the zero 
position, one pressure indicated by the dial was recorded 
along with that flow rate. The micromanometer was discon- 
nected from the pressure fitting and then adjusted to 
return the meniscus to its zero position. If the dial 
indication was again zero, the pressure difference obtained 
was considered to be correct. Occasionally the micromano- 
meter would not re-zero and the pressure measurement at 


that location was retaken after the instrument was adjusted. 
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The pressure drop over two foot portions of both the 
straight section and the curved section were taken to verify 
the other results. For measurement over these intervals 
the upstream pressure tap was connected to the top of the 
mMicromanometer. The downstream tap was connected to the 


reservoir of the micromanometer as in the previous setup. 


Pee RESSURE DROP 

The pressures measured were plotted downward on Figure 
31. The top of this graph represents atmospheric pressure, 
and the amount that the static pressure in the channel was 
below the atmospheric level can be seen. The pressures are 
plotted against pressure tap number. The 28 pressure taps 
were three inches apart along the channel, and numbers one 
mamough 16 were on the straight portion of the channel. 
The pressure for six different Reynolds numbers were plotted, 
and the same patterns can be seen in each. The irregulari- 
ties in the curves around pressure taps 14 and 15 and around 
taps 18 and 19 were of a smaller magnitude than the uncer- 
tainty for these measurements. Pressure taps 14 and 15 were 
mmenesstraight portion of the channel and no cause for the 
irregularities were apparent. However, burs or other dis- 
Eurbances which were undetected during construction of the 
channel may have caused these irregularities. Pressure taps 
eeand 19 were in the curved portion of the channel. In 
this case, the cause for the pressure irregularities may 
have been a difference in the position of the pressure taps 


relative to the vortices. A pressure tap that was being 
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affected by the stagnation of outward moving fluid and one 
under the influence of moving fluid would probably give dif- 
ferent pressure indications. Further study would be required 
to confirm that this was the cause of the irregularities. 

The slope of the pressure curve in the straight section 
should compare with the exact analytic solution for laminar 
flow between parallel flat plates. The results of such an 


analysis [Ref. 6] show that the pressure drop is 


where a is half the plate spacing and uw is the viscosity of 
the fluid. U is the mean velocity of the parabolic velocity 
profile between the parallel flat plates. However, the 
experimentally determined U was the volumetric flow rate 
through the rectangular channel divided by the cross sec- 
tional area. For a rectangular channel of aspect ratio 40, 
the difference in these two values of U would be small. 

The analytic value.of U for a rectangular channel could be 
obtained from an infinite series given in the literature 
[Ref. 19]. The first term of this sees would reduce U 
for the rectangular channel to 98.5315 percent of the value 
for an infinite parallel plate construction. The remaining 
terms of this series would rapidly become smaller. The 
experimental uncertainty was larger than this analytic 
difference, so that no correction was applied to the theo- 


retical pressure drop. 
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The experimental pressure drop was calculated by a com- 
puter program for both the straight and the curved portions. 
This program fit a least square straight line to nine con- 
secutive pressures representing a two foot section of the 
channel. The straight section of the channel was represented 
by pressure taps three through eleven inclusive. This 
avoided both the entrance region and the irregularities that 
were noted before. The curved portion was represented by 
pressure taps 19 through 27, where the vortices were most 
developed. The slope of each line was considered as the 
pressure drop in that portion of the channel. This program 
and its results are shown in Appendix C. These pressure 
drops in inches of water per foot of length, along with the 
analytic pressure drop calculated with the experimental 
values of U are shown in Table I. The percentage increase 
in the pressure drop for the curved section was calculated 
with respect to both the theoretical value and the measured 
value of pressure in the straight section. These results 
were also Repaired am Table: 1. The negative values shown 
were probably due to experimental uncertainties. 

Because eee significance might be obtained from a com- 
parison of pressure drops in the straight section with 
pressure drops in the curved section, where the vortices 
were fully developed, the same computer program was used to 
Sonsider one foot section of the channel. The straight 
portion was represented by pressure taps six through ten, 


and the curved section by the last one foot, pressure taps 
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mamemrough 28. The computer results for these calculations 
PeesshOwn in Appendix C. The resulting pressure drops and 
percentage increases are shown in Table II. These results 
showed no negative values for the increase, and had the same 
magnitude as the other results. 

To verify and compare the other results, direct measure- 
ments of the pressure drop in the two foot intervals used 
above were made. The error in this type of measurement could 
have a magnitude as high as 0.002 inches of water. Never- 
theless, the results, along with the resultinc percentage 
increases in the pressure drop, are presented in Table III. 
The values of the pressure drop measured in the straight 
section of the channel compared moderately well with the 
theoretical values. This comparison further supported the 
methods used to obtain the data. 

The percentage increase in the pressure drop from the 
straight section to the curved section was arrived at with 
three different measurements, and calculated with both the 
theoretical and the measured values in each case. The six 
different sets of percentage increases were plotted against 
the Reynolds numbers in Figure 32. The average increase for 
each Reynolds number is also identified in Figure 32. These 
results have considerable scatter due to the effects of 
experimental uncertainty. With an increase in the Reynolds 
number, there was an increased influence of the Taylor- 
Goertler vortices on the pressure drop. Even though the 
secondary flow was present, the increase in pressure drop in 


the curved section was small at the lower flow rates. At 
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o least squares 2 ft. (calculated) 
9 least squares 2 ft. (measured) 
4 least squares 1 ft. (calculated) 
4&4 least squares 1 ft. (measured) 


e Measured values 2 ft. (calculated) 


Percentage Increase in Pressure Drop 


i a Measured values 2 ft. (measured) 
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Figure 32. Percentage Increase in Pressure Drop 


vs. Reynolds Number. 
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higher Reynolds numbers, where the vortices became more 
active, the pressure drop in the curved section increased 
more than in the straight section. 

This data was nondimensionalized as the pressure drop 
coefficient, which was defined as 


Cp = —E-; 
4p U 


where p was the pressure measured at each pressure tap. 
There was a larger error bound on these values hecause the 
uncertainty in U was added to that in the pressure drop. 

The calculation of the pressure drop coefficient was done by 
hand on a desk calculator. A sample calculation for the 
lowest Reynolds number is given in Appendix C. The results 


were plotted in Figure 33. 
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VI. HEAT TRANSFER 


A. EXPERIMENTAL PROCEDURES 

In order to determine the effect of Taylor-Goertler 
vortices on heat transfer it was necessary to measure the 
heat transferred from a uniform temperature surface in one 
wall of the straight portion of the rectangular channel. 
The heat transfer coefficients of interest were defined by 


the equation 
Q= H * A*AT 


where Q was the heat convected, H was the heat transfer 
coefficient and A was the area of the heated surface. The 
temperature difference, AT, used depended on the particular 
heat transfer coefficient to be found. For the overall 
heat transfer coefficient, the log mean temperature differ- 
ence was used. The mean heat transfer coefficient was 
defined with the difference equal to the wall temperature 
minus the inlet bulk temperature of the fluid. The temper- 
ature difference used for the local eee transfer coefficient 
was the difference between the temperature of the individual 
heating element and the bulk temperature of the fluid enter- 
ing the control volume below that element. The differences 
in these definitions and the areas used will become clear 
during the discussion of the results. 

A heater plate with 16 separate resistance heating 


elements placed along the downstream length of the heater 
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was designed and built as described earlier. To measure the 
heat transferred from a uniform temperature surface it was 
necessary to control the Veltage te edel Of the to heating 
elements. The resistance and voltage at each heating strip 
was to be measured to determine the power generated. The 
temperatures were to be measured at two locations on each 
heating elements. 

After the heater was in place in the channel wall, the 
resistance of each heating elements was measured with a 
Rosemont Commutating Bridge, Model 920A. The high current, 
regulated D.C. power supply fed the power controller which 
governed the voltage across each heating element. This 
specially constructed equipment had never been previously 
used or tested. The output level of each of the 20 circuits 
was set as low as possible. The input voltage was adjusted 
until the power controller's reference voltage was between 
the maximum recommended value of 8.800 volts and 8.700 volts. 
To insure that the heater plate did not become too hot, flow 
was established in the channel. The voltage in each of the 
eee 16 power circuits was slowly increased, and the tem- 
perature along the plate was checked frequently. At any 
given power setting, four to five hours were required for 
the heated to come to equilibrium. As expected, in order 
to maintain the same temperature for each strip, the circuits 
near the leading edge of the plate required more power than 
the other circuits. If the variation of the power settings 


was not correct, the temperature distribution approached by 
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the plate was not uniform. When it was determined that the 
temperature would not be uniform, the power to the low 
temperature heating elements was increased and another 
settling period was allowed. Considerable time was required 
@ementain the first uniform temperature distribution. There- 
after, to obtain higher temperatures across the plate, each 
‘potentiometer's adjustment screw was turned equal amounts. 
As equilibrium was approached, small adjustments were made 
to remove temperature irregularities along the plate. In 
this manner, a uniform temperature could be achieved in a 
minimum time of about six hours. 

During continued testing of the equipment, it was found 
that the heater plate could sustain temperatures as high as 
200° F without damage. Higher operating temperatures were not 
attempted. During the early testing of the power controller 
the power circuit voltages were found to fluctuate even 
when the reference voltage was steady. Examination of the 
power circuit voltage showed a 60 cycle A.C. voltage super- 
imposed on the controlled D.C. voltage. The source of this 
A.C. voltage was the grounded side of the regulated power 
supply. Consequently, downstream from the protective 
diodes, a 0.05 uf capacitive connection was made across the 
inputs to the power controller. This connection removed the 
A.C. voltage and the fluctuation in the power levels. When 
the power circuits were loaded, the stability was found to 
be approximately +.004 volts over a 30 minute period. The 
power available was actually more than required for these 


tests. 
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After the equipment was tested, data collecting began. 
At the beginning of each data run, the flow rate was estab- 
lished and power was supplied to the heater. After two 
hours, an ice bath was made and the temperatures along the 
eee s were periodically checked. During the third hour, 
adjustments were made in the temperature and flow rate. The 
reference voltage was also checked at this time so that its 
Variation could be determined later. 

The temperatures measured and recorded included the bulk 
temperatures of the entering and exiting air, the tempera- 
tures at the 32 thermocouple locations along the heater, and 
the temperatures of the four locations on the outside of 
the heater plate. The thermocouples were connected to a 
digital voltmeter by a set of switches. The switching and 
recording of the thermocouple outputs were accomplished by 
hand. As equilibrium was approached, the temperatures were 
recorded at 30 minute intervals. When the temperatures had 
stopped changing within the limits of accuracy of this 
experiment, the values were recorded as data. The reference 
voltage and flow rate were rechecked to insure no change had 
taken place. Then the digital voltmeter was connected to 
the power controller and the voltage across each heating 
element was recorded. The digital ee 1 ene aes reconnected 
to the switches. The temperatures were recorded one addi- 
tional time, and the flow rate was again checked. The data 
was considered valid if no changes had taken place. The 


process of taking data took approximately 30 minutes. 
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Bee REOULTS 

In order to calculate heat transfer coefficients from 
the data collected, it was necessary to account for the heat 
messes. in the geometry of this experiment, there were 
several potential losses. The largest was expected to be 
the conduction through the back side of the heater plate. 
For this reason, the insulation, described earlier was placed 
over the heated area. To obtain an estimate of the conduc- 
tion losses through the plate, a finite difference model of 
the experimental situation was used. 

A computer program was used to obtain a numerical solu- 
tion of the transient and steady state temperature distribu- 
tions and the heat flows of the three dimensional model. 

The program used is entitled TRUMP and was developed at the 
Lawrence Radiation Laboratory oes 23). This program was 
implemented for the IBM 360 computer at the Naval Postgrad- 
uate School by C. Erbayram [Ref. 24]. Because of the length 
of this program and its output, only the steady state results 
of one run are included in Appendix D. 

In the computer program, one half of the symmetric 
arrangement of the heater plate, the insulation, and the 
surroundings were modeled by a total of 197 node points and 
606 thermal connections. The initial temperature of all 
the nodes, the heat generated at the nodes representing the 
Meating elements, and the heat transfer coefficients repre- 
senting convection to the surroundings, were input data. 


The results were the steady state temperature at, and the 


89 





heat flow from each of the nodes. The time to reach steady 
state was also given by the TRUMP program, and compared well 
with the six hours required in the actual experiment. 

After the TRUMP program was tested, the heat generation 
rates actually measured were input to the program. The heat 
transfer coefficients representing the convection to the air 
in the channel were initially estimated and then iterated 
until the wall temperature arrived at by the computer program 
was the same as that actually measured. In this way, the 
TRUMP program was used to model each experimental run. The 
heat flows from the nodes representing the heating elements 
were used to calculate losses in the reduction of the exper- 
imental data. This process will be described next. 

Temperatures and several heat transfer coefficients were 
computed from the data collected by the computer program in 
Appendix D. A table of the thermocouple calibration data 
was included in the computer program. All of the tempera- 
tures were calculated by interpolating the thermocouple 
readings in this ete: The local temperatures along the 
wall aa the average wall temperatures are shown in the 
computer results contained in Appendix D. 

The log mean temperature difference was calculated from 


the equation 


Tho - Tbi 
a 7 imennon! 
nN fw - Tbo 


where Tw was the average wall temperature, and Tbi and Tbo 


were the bulk temperatures of the air into and out of the 
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Meevewocction respectively. Using a table of fluid prop 
erties vs. temperatures, the fluid properties were evaluated 
at the average bulk temperature. The overall heat transfer 


Seeericient was then calculated from the equation 


* * OY = - 
Ho = 2 Cpe be /= hoa) 


Aw * Tm 


where V is the volumetric hiewerate, and “Aw is the areaot 
the heated surface. The mean heat transfer coefficient was 
also calculated in the program. With the mean temperature 


Pmerrerence defined as 
Tm = Tw - Thi 


the mean heat transfer coefficient was calculated from the 


equation 


* * Y = ; 
um = & cp V (Tbo AV oyaty 


Ag om 


The symbols are the same as used for the overall heat trans- 
fer coefficient. These two values appear in the computer 
Sucput. 

To calculate the local heat transfer coefficients, the 
power generated at each heating ee was calculated from 
the resistance of the strip and the voltage across it. Two 
values of heat flow obtained from the TRUMP program were 
also read into the computer for each heating element. Since 
one node in the TRUMP program represented two adjacent heat- 
ing elements, the heat flow values from the TRUMP program 
were divided in half to represent the losses from the indi- 


vidual heating elements. At the end elements, the conduction 
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losses along the channel wall were attributed entirely to 
the first and last heating element. The two heat flows for 
each heating element represented the conduction loss and 

the heat convection into the channel. A loss factor was 
calculated by the computer program by dividing the conduc- 
tion loss by the sum of the heat flows. The power generated 
in each strip was multiplied by one minus the loss factor 

to give the useful heat at eachelement. 


The radiation loss was accounted for by the equation 


4 


Or =e * o0 * As * (Tw = Thi’) 


where € was the emissivity of the nichrome taken as 0.7 and 
where As was the surface area of each heating element, 
0.014783 square feet. The temperatures used were converted 
to degrees Rankine, and the Stefan-Boltzmann Constant, 0, 


10 


was equal to 0.199967 x 10 BTU/min a °R. ‘The radiation 


loss was subtracted from the heat generated at eachelement 
less the heat conducted out of the plate to give the heat 
convected into the channel. 

The last loss accounted for was the convection from the 
heated air in the channel to the surroundings. Two thermal 
resistances were calculated for this type of loss. The 
first included terms for the forced convection in the chan- 
nel, the conduction through the plexiglas lower Wale canG 
the free convection from the horizontal outside surface of 


the plexiglas. The second thermal resistance included terms 


for the forced convection in the channel, the conduction 
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through the aluminum side walls, and the free convection 
from the vertical outside surfaces. The forced convection 
was represented by the mean neat transfer coefficient 
calculated by the program. The free convection was repre- 
sented by correlation formulas given in the text, Heat 
feamerer, by J. Holman [Ref. 20]. The convection loss from 


the air through the walls of the channel was calculated from 


Ao* (Tp = 70-0) 


Qc = R 


where Tb was the bulk temperature of air entering the con- 
trol volume below each heating element, R was the thermal 
resistance for each case, and A was the corresponding area. 
The total loss of this type was the sum of the two Qc's 
calculated with each of the two thermal resistances. In 
summary the heat convected into the channel from each heat- 
ing element, designated as p(j), was the power generated 
at each element minus the conduction losses through the 
plate and the radiation losses. The heat which raised the 
bulk temperature ae the air, was this convected heat, p(j), 
minus the convection losses from the air. 

The new bulk temperature of the fluid after passing a 
heating element was therefore calculated by the computer 


Beeeraing to 


Tbh(jt+l) = Th(j) + ip(j)-Oc()) | 


p* cp* Vv 
The local heat transfer coefficient was then easily calcu- 


lated from the equation 


2 Fe 








che = famers Ve SENG er Gal 
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where As was the area of one heating element. The distance 
Semcach heating element from the leading edge of the plate 
was also calculated and printed out along with the power 
convected, the bulk temperature of the fluid, and the local 
heat transfer coefficient. The average of these local heat 
transfer coefficients was also displayed on the computer 
outputs found in Appendix D. 

Five data runs were taken and reduced in the manner 
described above. The five runs cover four different Reynolds 
numbers and a moderate range of wall temperatures. The 
resulting local heat transfer coefficients were plotted 
against the heating element numbers in Figures 34 through 38. 
The center of the first heating element was 0.137 inches 
from the leading edge of the plate, and each following ele- 
ment was 0.274 inches farther downstream. 

The local heat transfer coefficients tend to be large 
at the beginning of the plate, decreasing more slowly with 
increasing distance along the heater plate. There is con- 
siderable scatter on these graphs due to the experimental 
uncertainty. 

The heat transfer coefficient at the second heating 
element appeared consistantly low. This was due primarily 
to the method of calculating loss factors, which did not 
charge this heating element for any of the end losses. How- 
ever, insufficient information about the heat flow was 


available to improve the technique. Other factors affecting 
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Figure 34. 


Heat Transfer Coefficient Tw = 125.: 
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the scatter included the variation of temperature along the 
wall. This variation was as high as +2.0° F and was probably 
due to small irregularities in the surface of the plate. 

Gee first nichrome strip on the fifth heating element was 
known to be indented, and its temperature was usually 1.5 

to 2.0 degrees higher than the surrounding nichrome ribbons. 

Typical error bounds are included in each graph. Dis- 
pite the scatter, the decrease in the local heat transfer 
coefficient is very much like that for a flat plate. How- 
ever, the heat transfer coefficients are consistantly higher 
than those for a flat plate. 

A brief attempt was made to find some result in the 
iMeeerature with which to compare the heat transfer coef- 
meements. The majority of the literature for developing 
heat transfer in rectangular channels considers both walls 
being heated. The Nusselt number of thermally fully 
developed flow between parallel flat plates with one wall 
at constant temperature was found to be 4.0 [Ref. 16]. For 
this experiment, 5 Nusselt number of 4.0, based on hydraulic 
diameter, corresponds to a heat transfer coefficient of 
0.0244 BTU/min ft? °F. For the lowest Reynolds number 
tested, the local heat transfer approaches this value. For 
the other Reynolds numbers, the local heat transter Coerri— 
cient is higher than the fully developed value. 

This data was not nondimensionalized because it was not 
compared to any other results. The intended use GiCthas 


data was for comparison with similar heated plate data from 
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the curved section of the channel. The procedures used to 
obtain these results were considered adequate for comparison 
to later work. However, several possible improvements were 
identified by this testing. Such improvements would include 
designing a system with fewer heating losses and developing 
a construction technique with would result in a smoother 


heated surface. 
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VIE. CONCLUSIONS 


The observations made during this investigation have led 
to several conclusions about the nature of Taylor-Goertler 
vortices in a curved rectangular channel. The vortices were 
Seeerved .to.develop with distance along the channel. The 
importance of this observation was that an analytic descrip- 
tion of the flow would have to account for this growth. The 
vortices were observed to increase in amplitude and complex- 
ity with increasing flow rate. Thus at a higher Reynolds 
number, less distance along the channel was required to 
obtain the same degree of development. Several identifiable 
patterns of the developing flow were observed. The patterns 
were useful in determining the degree of development of the 
vertices. 

Maes study of the velocity profiles in the channel indi- 
cated that Taylor-Goertler vortices are a complex laminar 
flow. Bm cheribte, at flow rates below that necessary to 
Sustain turbulent flow, momentary bursting to turbulent 
flow was observed. 

The effect of Taylor-Goertler vortices on the pressure 
drop in the curved rectangular channel was observed to vary 
from no effect to approximately a 20 percent increase. The 
amount of the increase in pressure drop due to the Taylor- 
Goertler vortices increased with increasing Reynolds numbers. 


This seems reasonable considering that as the secondary flow 
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/mmereases in activity the interaction of the Reynolds 
stresses with the mean flow intensifies. 

The investigation of heat transfer from a constant 
temperature plate in one wall of the channel resulted in a 
procedure that would be adequate for determining the effects 
of Taylor-Goertler vortices on heat transfer. The heat 
transfer coefficients obtained for the straight section of 
the channel showed the trends expected and compared with 
the limiting data available. Improvements in the heat 


transfer measuring procedures could be made. 
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VIII. RECOMMENDATIONS 


At the end of this investigation, it was clear that 
Several recommendations could be made regarding continued 
Seveyeyor Taylor—Goertler vortices in a curved rectangular 
Channel. More flow visualization should be done to clarify 
the patterns in the vortex flow. Different aerosol injec- 
Piomecechniques, including a wall slot, should be tried. 
More intense lighting and perhaps narrower slits should be 
used to obtain clearer pictures of the flow. An improved 
Flow visualization technique should then be used to examine 
more locations along the curved portion of the channel. Due 
to the fact that there were motions in the flow patterns, 
still photography could not indicate all of the observations. 
A motion picture study of flow patterns would be ideal. 

More could be learned from a hot wire anemometer study 
where the probe could be moved sideways across the channel 
as well as Pea iainy in and out. In this manner, the perio- 
dic variation in the spanwise direction could be determined. 
The size of the vortices and their locations could then be 
determined. In addition, by knowing where the probe was 
in a vortex, the magnitude of the other velocity components 
could be determined. 

With respect to pressure drop, the placement of more 
pressure taps along the last portion of the straight section 


and along the entire curved section would improve the 
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reliability of the average pressure drop obtained. An 
investigation should be made to determine if the pressure 
variation in the direction perpendicular to the flow could 
be measured. 

It was clear that work should continue towards deter- 
mining the effect of Taylor-Goertler vortices on heat trans- 
beGeesAwheater plate for the curved section could be built 
in the same manner as the flat heater plate used in this 
experiment. However, it is recommended that a casting 
material with less shrinkage and better thermal insulation 
properties be used. If this is done, a smoother heated 
surface and fewer thermal losses would result. It is 
recommended that the development of a good heat transfer 
measuring technique be given considerable thought. 

For any analytic solution attempted, it was felt that 
it should account for the development of the vortices with 
distance and flow rate. If the magnitudes of the velocity 
components are determined, as suggested above, they should 
compare with any ee solution obtained. The inter- 
action of the Reynolds stresses with the mean flow would 
have to be accounted for to predict the increase in pres- 
sure drop. Only a non-linear analysis such as J. Stuart's 
work [Ref. 25] could be expected to give such results. 
Other analytic techniques including numerical methods ought 


Gemme tried. 
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APPENDIX A 


ERROR ANALYSIS 


A calculation of the uncertainty for each of the major 
variables of this experiment was made in accordance with 
the methods of S. Kline and F. McClintoch [Ref. 26]. The 
estimates of the uncertainty in the measured quantities were 
made quite conservative so that there was considerable con- 
fidence in the uncertainties arrived at. As an example of 
the calculations the uncertainty in Reynolds number is 


calculated below. The Reynolds number is given by 


Ss DG 
v 


Re = 
and the uncertainty is taken as 


dre _ / au“, @ph*_ av” 


= ——s <a + ee 


Re U Dh Vv 


The uncertainty in the flow velocity is obtained from other 
calculations like this one based on estimates of the uncer- 
tainty in reading the volumetric flow rate through the roto- 
meter and the uncertainty in the cross sectional area of the 
channel. The results of those calculations were oF = 0.02805. 
The uncertainty in hydraulic diameter and kinematic visco- 
Sity are obtained from estimates as 0.00426 and 0.00136 
respectively. Then the uncertainty in the Reynolds number 


is 
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dR = 
cee = \/ (7.868 + 0.1815 + 0.0185) x 10°" = 0.02840 


Other quantities and their uncertainties are given below. 
Mimemerror bound for any quantity can be obtained by multi- 


plying its magnitude by its uncertainty. 


Quantity Uncertainty 
Ac 0.00412 
Aw ' 0700320 
Cp 0. 07sl5 
cp 0.00010 
De 0.0 Zee 
dH 0.00426 
H O22 750 
Ho 0203 st 
p (pressure) 0.05000 
Ph (heat generated) O20 SoZ 
p(j) (heat convected) Onz5502 
Re | 0.02840 
T 0.00750 
T&m C-0 2000 
Tw On0 260 7 
U 0.02805 
V O02 775 
fe) 0.00080 


107 





APEENDLX. B 


COMPUTER RESULTS FOR VELOCITY AND TURBULENCE 


The computer program for reducing the hot wire anemometer 
calibration data is shown below. This program was adopted 
from a program written by Mr. Loyd Smith while he was study- 
ing at the Naval Postgraduate School. The resulting cali- 
bration curve is also shown. 

The hot wire anemometer data reduction program was also 
adapted from a program by Mr. Smith and is shown below. The 
reduced data for the seven velocity and turbulence profiles 


Pomenen included. 
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APPENDIX C 


CALCULATIONS AND COMPUTER RESULTS FOR PRESSURE DROP 


The computer program written to find the slope of the 
pressure curve is shown below. The results for six dif- 
ferent Reynolds numbers are shown. The first set of results 
represent nine pressure readings or two foot intervals. The 
second complete set represents selected one foot intervals. 
The last item shown in this appendix is a tabular calcula- 
tion of the pressure drop coefficient for the lowest Rey- 


nolds number tested. 
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SLOPE OF PRESSURE DROP CURVE FGR REYACLES NUMBER = 609.0 


POSITION (FT) FRESSURE (IN F2C) 
0.0 Ox 
0.250GC GlOL4 so 
0.50000 C.CLSCGC 
0.750CC 0.01650 
1.00000 0.018C0 
THE PRESSURE DROP IN THE STRIGHT SECTICN IS 0.004060 
SLCPE OF PRESSURE DROP CURVE FOR REYNCLES NUMBER = 606.0 
POSITION (FT) PRESSURE (IN H2C) 
0.0 Ceeagac 
0.25000 C1 0aGeC 
C.500CC C.04000 
0.75006 C.04150 
1.000CO0 0.0420C 
THE PRESSURE DROP IN THE CURVED SECTION IS 6. 0CsCec 
SLCPE GF PRESSURE DROP CURVE FOR REYACLCS NUMBER = 875.0 
POSITION (FT) PRESSURE (IN 20) 
0.0 C.01950 
C.250C0 0.02100 
0.5000C Cece 2 G0 
0.75000 0.025CC 
1.000CGC 0.02700 
THE PRESSURE DRGP IN THE STRIGHT SECTION IS 0.007600 
SLOPE OF PRESSURE DROP CURVE FOR REYNCLES NUMBER = 875.0 
POSITION (FT) FRESSURE (IN F2C) 
One 0.05000 
0.25000 C.052CC 
0.50000 C6 C54CC 
0.7500C : 
1.00000 C.058C0 
THE PRESSURE DROP IN THE CURVED SECTICN IS 0.0080C0 
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SLOPE OF PRESSURE DROP CURVE FOR REYNCLES NUMBER = 1166.0 


PCSITION (FT) PRESSURE (IN H20) 
Q.0 . 
0: 25000 5203300 
Q0.500CC Q0.C3500 
Q0./5000 C02 UCC 
1.0Q000CC QO .03900 
THE PRESSURE DROP IN THE STRIGHT SECTION IS C.CC84CC 


SLOPE OF PRESSURE DROP CURVE FOR REYACLES NUMEER = 1166.0 


POSITION (FT) PRESSURE (Clie 2u)) 
0.0 0.07100 
0.25000 C.0/73C0 
0.50000 C.C/éECC 
O. 75000 0.0/900 
L.000CQ O.C82C0 
fee PRESSURE DROP IN THE CURVED SECITTCN Es O.01L1L2C0 


Seeerc OF PRESSURE DROP CURVE FOR REYNELCS NUMBER = 1484.0 


POSITION (FT) PRESSURE (CIN H2aGy 
0.0 0.04000 
0.25000 C.043C0 
0.50000 Q0.C4E0C 
O.750CC 0.04SC0 
1.00000 C.C52CC 
Weer RESSURE DROP IN THE STRIGHT SECTIEGNETS C.0120C0 


SLCPE CF PRESSURE CROP CURVE FOR REYACLCS NUMBER = 1484.0 


POSITION (FT) PRESSURE (CIN h20) 
Q0.0 C.095CC 
0.25000 C.098CC 
O.500CC O.10200 
QO./5000 C.1LO5CC 
1.00000 0.10900 
THE PRESSURE CROP IN THE CURVED SECTION IS 0.014CCO0 
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SLOPE OF PRESSURE OROP CURVE FOR REYACLOS ALMEER = 1748.C 


POSITION (FT) PRESSURE @CiiN Ee) 
0.0 ° 
0% 25000 620320 
° C.O057TC 
Oe oCOC 6108100 
1.00000 C.C&4CC 
Mite PRESSURE DROP IN THE STRIGHT SECTICN IS 0.0144C0 


SLCPE OF PRESSURE DROP CURVE FOR REYNCLCS NUMBER = 1748.0 
POSITION (FT) PRESSURE (IN H20) 
Ono. Gz VIZCo 
0.25000 0.116CC 
0.500CC 0.2100 
6.75000 C.12ECC 
1.00000 0.13C0C 
THE PRESSURE CROP IN THE CURVED SECTION IS C.C184CC 
SLCPE OF PRESSURE DRCP CURVE FOR REYACLCS NUMBER = 1987.0 
POSITION (FT) FRESSURE (IN #20) 
0.0 C.06CCC 
0.250CC C.06400 
0.5000C C.068CC 
0.75000 0.07206 
1.000CC 0.07600 
THE PRESSURE DRGP IN THE STRIGHT SECTICN IS C.01€4CC 


SLOPE GF PRESSURE DROP CURVE FOR REYNCLOS NUMEER = 1987.0 


POSITION (FT) PRESSURE MCN i2u 
0.0 0.13200 
0.25000 C.137C0 
0.50000 C.142C0 
QO. 75000 0.14/7C0 
1.000C0 C.152C0 
THE PRESSURE DROP IN THE CURVED SECTICN IS 0.0200C0 
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APPENDIX D 


COMPUTER RESULTS FOR HEAT TRANSFER 


The final results for one run of the TRUMP program are 
shown on pages 130 through 144. The computer program for 
the calculation of the heat transfer coefficients follows 
the TRUMP program. The five sets of heating data results 
are presented with the temperature distribution shown on 
the first page and the heat transfer coefficients on the 


second page of each of the results. 
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